Scaling and Decoherence in the Out-of— Equilibrium Kondo Model 
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We study the Kondo effect in quantum dots in an out-of-equiiibrium state due to an appiied 
dc-voitage bias. Using the method of infinitesimal unitary transformations ("flow equations"), 
we develop a perturbative scaling picture that naturally contains both equilibrium coherent and 
non-equilibrium decoherence effects. This framework allows one to study the competition between 
Kondo effect and current-induced decoherence, and it establishes a large regime dominated by 
single-channel Kondo physics for asymmetrically coupled quantum dots. 



Since the first experimental observations of the Kondo 
effect in the Coulomb blockade regime of quantum dots 
IH ) a wealth of experimental and theoretical work 
has addressed the properties of this highly controllable 
correlated electron system. If a quantum dot weakly cou- 
pled to two leads carries a net spin, resonant tunneling 
through the dot becomes possible and leads to a Kon- 
doesque increase of the conductance up to the unitarity 
limit upon lowering temperature 0, S 0- These real- 
izations of the Kondo effect in quantum dots have led 
to new questions related to the out-of-equilibrium na- 
ture of the Kondo system with a stationary current for 
an applied voltage bias. Despite many theoretical efforts 
(e.g., Refs. BIlMHmUlEi), a satisfactory theory 
for the out-of-equilibrium Kondo effect does not yet ex- 
ist. Most theoretical methods that have been developed 
for equilibrium [J^ cannot easily be generalized to the 
non-equilibrium situation. 

In this Letter, we focus on the case of large voltage 
bias V Tk at zero temperature (T = 0) as a step 
towards a more complete understanding of the out-of- 
equilibrium Kondo model. Kaminski et al. first 
suggested a "poor man's" scaling method, and subse- 
quently Rosch et al. 01 developed a more sophisti- 
cated approach based on frequency-dependent vertices 
and Keldysh diagrammatics. Both groups noted that de- 
coherence generated by the current is essential since it 
introduces the decoherence rate r rc i oc V/ In (V/Tr) due 
to non-equilibrium spin relaxation processes. However, 
a scaling picture in a Hamiltonian framework in which 
this decoherence scale emerges naturally has as yet not 
been developed. This Letter shows that the method of 
infinitesimal unitary transformations ("flow equations") 
[lil Hif provides such a suitable generalization of Ander- 
son's "poor man's" scaling picture 17]. It allows one 
to study the flow of coupling constants and the ensuing 
phase diagram in a way that is similar to the scaling anal- 
ysis of an equilibrium problem, and it establishes a large 
regime dominated by single-channel Kondo physics for 
asymmetrically coupled quantum dots. 

We consider a Hamiltonian describing a spin-1/2 de- 
gree of freedom S coupled to two leads a, a 1 = l,r with 



voltage bias V as an effective model for a quantum dot 
in the Kondo regime 
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Here s (a , p , )(ap) = \ Y, a ,0 c 

tum labels and fj,^ r — ±V/2. If the quantum dot can be 
described by an Anderson impurity model with tunneling 
rates T[ r from the left /right lead, the coupling constants 



are related by J 2 



and Ju/Jrr = L;/F r (notice 



J lr — J rl for hermiticity) fic| . 

The flow equation method makes a Hamiltonian 
increasingly diagonal by applying a sequence of infinites- 
imal unitary transformations with an antihermitean gen- 
erator rj: 4§ = [r)(B),H(B)]. Here B labels a one- 
parameter family of unitarily equivalent Hamiltonians 
and has the dimension (Energy) -2 : H(B = 0) is the 
initial Hamiltonian and H{B) is the unitarily equiva- 
lent Hamiltonian with matrix elements with energy dif- 
ferences \AE\ > _B -1 / 2 being eliminated. This RG- 
like separation of energy scales for B — > oo can be 
achieved by choosing the generator as the commutator 
of the diagonal and the interaction part of the Hamil- 
tonian: rj(B) = [Hq(B), H- m t(B)]. In contrast to con- 
ventional scaling approaches that successively eliminate 
high-energy states in the Hilbert space, the flow equation 
method keeps all states but makes the scattering pro- 
cesses increasingly more energy-diagonal. This method 
has been successfully applied to numerous equilibrium 
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FIG. 1: Left: Conventional scaling picture where states are 
integrated out (here depicted for cutoff A smaller than the 
voltage bias). Right: Flow equation approach. Here all scat- 
tering processes with energy transfer |A£| < _B -1 / 2 are re- 
tained in H(B). 
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many-body problems (e.g., Refs. 0,0, Hjj), where it 
correctly describes the RG-flow in which the UV-cutoff 
is identified as A oc B^ 1 / 2 . For an out-of-equilibrium 
model like the Kondo model with voltage bias, the differ- 
ence between state elimination and flow equation diago- 
nalization turns out to be more fundamental (see Fig. [T|> : 
in the flow equation picture scaling below B~ x l 2 < V is 
straightforward. In particular, the flow equation Hamil- 
tonian H(B) for B~ x l 2 < V still describes the stationary 
current flowing across the dot since energy-diagonal scat- 
tering processes are not eliminated, which is essential for 
obtaining the current-induced decoherence scale. 

We apply the flow equation approach with the canon- 
ical choice of T) (where Hq is the kinetic energy of 
the electrons) to the Kondo Hamiltonian Q. During 
the flow higher order interactions are generated and we 
parametrize the Hamiltonian as 



H(B) = J2 e t4 a c ta + Y, J t't( B )S-s t , t 



(2) 
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where : . . . : denotes normal-ordering with respect to 
the unperturbed Fermi sea fl5| . Here t' , t are general in- 
dices, initially K vt U , U (B = 0) = and we neglect newly 
generated normal-ordered terms in 0(J 3 ) and higher. 
Lengthy but straightforward calculations lead to the fol- 
lowing set of flow equations |2*l|: 
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Here the Fermi sea expectation values n + (it) 



dof , 4. 
= ( r I / 



and n~(u) A = (c u cl) arise due to the normal-ordering 
prescription. A nontrivial test for this calculation is 
provided by the equilibrium case where t just describes 
the momentum label. Using the approximate (however, 
asymptotically correct in the IR-limit) parametrization 
pJp>p(B) = g(B) exp(— B(e p ' — e p ) 2 ) with the dimension- 
less IR-coupling constant g(B) = pJ tF<iF {B) in Eqs. 

one derives the conventional third order scaling 
equation for the equilibrium Kondo model, dg/dhiA = 
-9 1 + .9 3 / 2 (with the identification In A = -(1/2) \\\B). 

Eqs. ©, I0J contain complete information about the 
Hamiltonian flow with voltage bias to 0(J 4 ),0(J 3 ), 



resp., and can be analyzed without further approx- 
imations |2l| . However, in order to gain analyt- 
ical insight one can employ the following approxi- 
mate parametrization that focuses on the IR-limit: 
pJ{ap')(a P ){B) = g a (B)exp(-B(e P ' - e p ) 2 ) with g a (B) = 
pJ(ap, a )( al j, a )( B ) for a = l,r, and pJ ( i p r ){rp) {B) = 



pJ(r P )(ip')(B) = g t (B)exp(-B(e p > 



with gt(B) = 



l/(pi - p r ) fj£ dtpJ( le )( re )(B). Here gi and g r are the 
coupling constants for left-left and right-right scattering 
processes located at the respective Fermi surfaces. For g t 
we choose an average over the transport couplings since 
this average is directly related to current and conduc- 
tance (2^. Inserting these parametrizations into Eqs. 
(@J , one arrives at the following set of equations that have 
to be integrated starting from B = D~ 2 (D is the initial 
UV-band cutoff, a = l,r): 
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where dkb c /dB = gbg c with the initial condition kb c (B = 
D~ 2 ) — for all b,c = I, r, t. The scaling picture deduced 
from Eqs. (JSJ, J^J is the main result of this Letter. 

Eqs. 0, JBJ take different forms for B^ 1 / 2 > V and 
B~ x l 2 < V. We first analyze the scaling behavior down 
to the scale set by the voltage bias; only terms up to 
second order need to be taken into account here (higher 
order terms are unimportant for V 3> Tk). One arrives at 
a set of equations already analyzed in Ref. 10] (a = I, r): 



dga 
dB 



2B 



(9 2 a+9t) 



dgt 
dB 



2B 



9t(9l+9r) ■ (7) 



In the following discussion we only look at Kondo dots 
described by an Anderson impurity model. The scaling 
invariant Kondo temperature is then set by the equilib- 
rium case, Tk = D^fgi + g r exp(—l/gi+g r ), and one eas- 
ily shows g t (B = V- 2 ) = (VWr/Ti + T r ) / ln(V/T K ). 
From © one deduces that the growth of gt effectively 
stops below the voltage bias scale since there is no sharp 
Fermi surface for transport processes. However, looking 
at the second order terms in (0 the strong-coupling di- 
vergences for gi and g r are not cut off. This has led to the 
prediction of two-channel Kondo physics in Kondo dots 
with voltage bias (23|. But even in the weak-coupling 
regime third order terms in the coupling constant eventu- 
ally become more important than second order terms for 
nonvanishing voltage bias in (J5J, @. For B~ x l 2 < V the 
dominant terms in the flow equations are approximately 
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(but sufficiently accurate for a qualitative picture) 

dB-~7^y8 9t (8) 

dg a 1 / 2 , o d\ng t \ , , fn s 

-dB = 2B{ 9a + 29a d^B) fOTa -^- W 

The flow changes qualitatively below the decoherence 
scale B" 1 / 2 = Vfl£(B = V"- 2 ): for B > B dcc algebraic 
decay gt(B) oc i? -1 / 4 sets in. In Eq. 10 one can then 
study the competition between coherent strong-coupling 
behavior from the second order term, and decoherence 
effects that arise in linear order in g a for B > i?d cc - The 
growth of the coupling constants gi , g r stops at the deco- 
herence scale unless the coupling constants have already 
become too large. This qualitative analysis is confirmed 
by the numerical solution of the full set of equations 
©, © depicted in Fig. |2 Also notice that current- 
induced decoherence enters differently from temperature 
into the flow equation (JHJ. Temperature acts as an in- 
frared cutoff in the Kondo strong-coupling terms g^/2B, 
whereas current -induced decoherence and the coherent 
strong-coupling processes are in competition. Therefore 
the Hamiltonian flow derived in the Letter confirms the 

— 1/2 

existence of the decoherence rate IY^i oc -B dec due to 
current-induced spin relaxation 0, [n| . It is essential 
to work in a framework where energy-diagonal processes 
are retained because this means that a "window" of order 
voltage bias is open for transport (compare Fig.^l. It is 
exactly these transport processes that are responsible for 
the emergence of the decoherence scale since they lead to 
the V r i? _1 / 2 -terms in (JHJ and 0. A conventional scal- 
ing approach that removes states around the two Fermi 
surfaces and therefore purports to treat energy-diagonal 
transport processes (see Fig. [TJ cannot describe this and 
leads to 1/B-terms instead j24j . 

Since we now know explicitly how decoherence enters 
into flow equations J^J, ©, we can determine a quanti- 
tative phase diagram. The scaling curves for symmetri- 
cally coupled Kondo dots in Fig. |2Ji confirm the absence 
of two-channel Kondo physics for large voltage bias in 
the sense that all couplings remain small. We define 
the strong-coupling regime by requiring that at least one 
coupling grows larger than 0.5; the resulting curve is de- 
picted in Fig. |21 25] . This line should not be interpreted 
as a phase transition; the crossover is expected to be 
smooth — similar to the effect of temperature in the equi- 
librium Kondo model. For asymetrically coupled Kondo 
dots this strong-coupling regime extends to remarkably 
large values of the voltage bias since the decoherence 
scale -B dcc ' is proportional to the current I (see below) ; 
for a given ratio of V/Tk the maximum value of I/Tk 
is achieved for symmetrically coupled dots (see below). 
Therefore current-induced decoherence is less effective 
in asymmetrically coupled Kondo dots when competing 
with the coherent inter-lead strong-coupling processes. 
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FIG. 2: Universal curves for the flow of the coupling constants 
gi (full lines), gt (dashed lines) and g r (dotted lines) for sym- 
metrically coupled quantum dots (Plot a with Fi/F r — 1) 
and for an asymmetrically coupled quantum dot (Plot b with 
Ti/Tr = 4). Results are shown for various ratios V/Tk la- 
belling the resp. curves from top to bottom, gi and g r coincide 
in the symmetric case, therefore only gi is shown in a. 



Experimentally, this phase diagram can be explored by 
measuring the Kondo dot density of states Pd(^), e.g. via 
a 3-lead setup |2rJ ] . The strong-coupling regime in asym- 
metric Kondo dots implies a density of states at the Fermi 
level of the more strongly coupled lead that only drops 
significantly (< 25%) below the Friedel value once V/Tk 
is well in the weak-coupling regime in Fig. [3] For already 
small asymmetries the crossover to the strong-coupling 
regime is driven by single channel Kondo physics in the 
sense that the couplings at one Fermi surface start to 
dominate strongly. This effect can be traced into the 
weak-coupling regime by the observation that the ratio of 
the local density of states at the two Fermi levels is then 
given by Pd(ftiX/Pd(^r) = " 2 with a = gi(Bi)T r / g r (B r )Ti 
from Fig. 01,27|. This observation Pd(^i) / Pd(p-r) 7^ 1 in 
asymmetric Kondo dots has also been made in other ap- 
proaches though without obtaining a quantitative phase 
diagram (compare Refs. flU Eil |28| h 

A final remark regarding the calculation of observ- 
ables: these need to be unitarily transformed as well 
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FIG. 3: Phase diagram of the non-equilibrium Kondo model 
as a function of asymmetry Ti /Y r and voltage bias. The full 
line separates a weak-coupling regime (defined by all cou- 
plings remaining smaller than 0.5 during the entire flow) from 
a strong-coupling regime. The contour lines (dashed) show 
the ratio a = gi(Bi)T r / g r (B r )Ti, see text (B a is defined as 
the £?-value where g a (B a ) takes its maximum value). 

(lif (they typically change their form completely once 
B > -Bdcc)- For example the current operator remains 
form-invariant with the flowing coupling gt(B) up to 

scale B doc : 1(B) = ig t (B)J2 p > <p S ■ (s*(; p ')(r P ) ~ S(r P >)(i P ))- 
One can then work with the renormalized effective Hamil- 
tonian on the scale -Bdoc and a conventional Keldysh 
calculation yields I = (37r/4)V~g 2 (-Bdcc)- This leads 
to the well-known perturbative result for the conduc- 
tance 10] G(V) = G u (37r 2 /16)/ln 2 (V r /T K ) where G u = 
(e 2 /'kK) 4r;r r /(r/-|-r r ) 2 is the conductance in the unitar- 
ity limit. Notice that transport quantities are not low- 
energy properties like the Kondo dot density of states 
at the Fermi levels and are therefore unaffected by the 
strong-coupling physics in Fig. as long as V 3> 7k- 

Summing up, we have developed a Hamiltonian scal- 
ing picture for Kondo dots with voltage bias that allows 
us to express physical quantities in terms of renormal- 
ized parameters. We confirm the absence of two-channel 
Kondo physics for symmetrically coupled quantum dots, 
and show the existence of a large regime dominated by 
single-channel strong-coupling physics for asymmetric 
dots. The flow equation approach presented here should 
be useful for other non-equilibrium models as well. 
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